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Recently much attention is paid to the role of the orbital degrees of freedom in transition metal 
oxides as it remains unclear whether they can remain in a quantum disordered state at zero tem- 
perature. Discrete symmetry of the orbital sector counteracts the quantum melting, but especially 
in doped systems there are signs of dynamical frustration involving the spin-, charge-, and orbital 
sector simultaneously. We discovered that even the simple Kugel-Khomskii model, describing e 9 
degenerate Mott-insulators, is characterized by a point of perfect dynamical frustration on the clas- 
sical level, reached in the absence of Hund's rule and electron-phonon couplings. This frustration 
is lifted on the quantum level, and the true nature of the ground state is still unknown. At present 
there are two proposals: the KCUF3 phase, stabilized by an order-out-of-disorder mechanism, or 
spin-orbital valence bond phases. It will be argued that at least in the Cu based systems of this 
kind, the electron-phonon coupling is primarely responsible for driving the systems away from the 
special point in the phase diagram. 
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I. INTRODUCTION 

The paradigm of correlated electron physics is based on 
the idea that for a certain category of systems one bet- 
ter starts out with the electronic structure of the atoms, 
treating the delocalization of the electrons in the solid as 
a perturbation. Any student of physics has to struggle 
through the theory of atomic multiplets, which is rather 
complicated because of the intricacies associated with 
orbital angular momentum. At first sight it is there- 
fore remarkable that these orbital degrees of freedom are 
completely neglected in the main stream of correlated 
electron physics. Recently the interest in 'orbitals' has 
been reviving, especially since they appear to be relevant 
in one way or another in the colossal magnetoresistance 
(CMR) manganites. In the wake of this development, 
questions are asked on the relevancy of these orbitals in 
the context of seemingly settled problems like the metal- 
insulator transition in V2O3 Qj. In this contribution we 
will review yet another recent development. Even in the 
Mott-insulating limit, where the physics simplifies con- 
siderably, the interplay of orbital and spin degrees of 
freedom poses a problem of principle. 

There are two limits where the role of orbital degen- 
eracy is well understood: (i) The 'band structure limit', 
which is based on the assertion that electron correlations 
can be neglected. In any modern local density approx- 
imation (LDA) band structure calculation, orbitals are 



fully taken into account on the one particle level, in so 
far as the atomic limit is of any relevance. These trans- 
late into various bands, giving rise to multi-sheeted fermi 
surfaces, etcetera, (ii) The localized, orbital and spin or- 
dered case which we will refer to as the 'classical limit'. In 
Mott-insulators, orbital degrees of freedom acquire a sep- 
arate existence in much the same way as the spins of the 
electrons do. The orbitals can be parametrized by pseu- 
dospins and these form together with the physical spins a 
low energy sector which is described by generalizations of 
the Heisenberg spin-Hamiltonian 0j . These are the spin- 
orbital models, like the Kugel-Khomskii (KK) model for 
e g degenerate cubic cuprates ||. The 'classically' ordered 
states, becoming exact in the limit of infinite dimensions 
(d — > 00) and/or large (pseudo) spin (S — > 00), define 
what is usually meant with orbital and spin order. 

The question arises if there are yet other possibilities. 
We started to study this problem quite some time ago 
Q, well before the subject revived due to the mangan- 
ites. Our motivation was actually related to a theoreti- 
cal development flourishing in the 1980's: large N theo- 
ries By enlarging the symmetry, say from SU(2) to 
SU(N) with N large, new saddle points (ordered states) 
appear which correspond to the fluctuation dominated 
(non-perturbative) limit of the large S /large d theories. 
For a single correlated impurity, orbital degeneracy leads 
in a natural way to these large N notions. We asked the 
question if these large TV notions could become of rel- 
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evance in lattice problems. We focussed on the simple 
problem of the e g Jahn- Teller degenerate Mott-insulator, 
rediscovering the KK Hamiltonian ||] . We tried to tackle 
this problem using the techniques invented by Arovas and 
Auerbach for the SU(N) symmetric Heisenberg model 
|@|. We found that the SU(4) symmetry is so badly bro- 
ken that the large N techniques were of little help, which 
is another way of saying that the physics of the KK model 
is not controlled by large global symmetry. However, we 
did find a special approximate solution which revealed 
that the quantum fluctuations are actually enhanced, and 
this motivated us to study these fluctuations in more de- 
tail starting from the large S limit. In this process we 
discovered that the enhancement of the fluctuations is 
due to the control exerted by a point in parameter space 
which can be either called an infinite order quantum- 
critical point, or a point of perfect dynamical frustration 
in the classical limit 0. 

This phenomenon will be discussed in the next sec- 
tion. It poses a rather interesting theoretical problem. 
So much is clear that the ground state degeneracy of 
the classical limit is lifted by quantum fluctuations and 
the question is on the character of the true ground 
state. As will be discussed, either the classical spin- 
orbital order might survive, stabilized by an order-out-of- 
disorder mechanism, or quantum-incompressible valence- 
bond like states might emerge. In Section III the role of 
electron-phonon coupling will be addressed, emphasizing 
the rather counter-intuitive result of LDA+U electronic 
structure calculations that phonons play a rather sec- 
ondary role despite the fact that the lattice deformations 
are large. Finally, the situation in the manganites will 
be shortly discussed in Section IV. 



II. THE KUGEL-KHOMSKII MODEL AND 
DYNAMICAL FRUSTRATION 



| , \z) = I l J . Pauli matrices a u (u = x, y, z) 
are introduced acting on these states. Different from the 
spins, the SU(2) symmetry associated with the pseu- 
dospins is badly broken because the orbitals communi- 
cate with the underlying lattice. Although the e g states 
are degenerate on a single site, this degeneracy is bro- 
ken by the virtual charge fluctuations, which take place 
along the interatomic bonds, i.e., in a definite direction 
with respect to the orientation of the orbitals. It is there- 
fore convenient to introduce operators which correspond 
to orbitals directed either along or perpendicular to the 
three cubic axes a = a, b, c, given by (tJ*— |) and (t?+^), 



where 



a(b) 
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In terms of these operators, the Kugel-Khomskii Hamil- 
tonian can be written as (J = t 2 /U and t is the hopping 
along the c-axis) , 
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neglecting the Hund's rule splittings cx Jjj of the inter- 
mediate d s states (Jh is the singlet-triplet splitting). In- 
cluding those up to order rj = Jh/U yields in addition, 



H 2 = JvJ2 
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Consider a Mott-insulator which is characterized by or- 
bital degeneracy, besides the usual spin degeneracy. Dif- 
ferent from pure spin problems, these spin-orbital prob- 
lems are rather ungeneric and depend on the precise sys- 
tem under consideration. A simple problem is a cubic lat- 
tice of 3(i-ions in a d 9 configuration: the Kugel-Khomskii 
problem, which directly applies to Cu perovskites like 
KCuF 3 or K 2 CuF 4 §]. The large Mott gap in the charge 
excitation spectrum simplifies matters considerably and 
one derives an effective Hamiltonian by insisting on one 
hole per unit cell, deriving superexchange-like couplings 
between the spin and orbital degrees of freedom by inte- 
grating out virtual charge fluctuations. 

The spins are described as usually in terms of an su(2) 
algebra (Si). The orbital degrees of freedom are the 
e g cubic harmonics x 2 — y 2 ~ \x) and 3z 2 — r 2 ~ \z), 
which can be parametrized in terms of pseudospins as 



Eq.'s are rather unfamiliar: they describe a regu- 

lar Heisenberg spin problem coupled into a Pott's like 
orbital problem (choose two out of three possibilities 
~ x 2 — y 2 , ~ y 2 — z 2 , ~ z 2 — x 2 ). 

The oddity of Eq. 's (0,||) becomes clear when one stud- 
ies the classical limit. As usually, the S"s and the r's are 
treated as classical vectors. In order to draw a phase 
diagram we introduced another control parameter, 

H 3 = -E z J2t*, (4) 



a "magnetic field" for the orbital pseudo-spins, loosely 
associated with a uniaxial pressure along the c-axis. The 
classical limit phase diagram as function of rj and E z is 
shown in Fig. 1. 
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FIG. 1. Phase diagram of the Kugel-Khomskii model in 
the classical limit, as function of the Hund's rule coupling Jh 
and tetragonal crystal field E z (reproduced from Ref. |7J). 

For a detailed discussion of the various phases we refer 
to Ref. [Q. To give some feeling, for large positive E z the 
x 2 — y 2 orbitals are occupied, forming (a, b) planes of anti- 
ferromagnetically coupled spins (AFxx). This is nothing 
else than the situation realized in, e.g. La2Cu04. For 
large negative E z the 3z 2 — r 2 orbitals condense, forming 
a 3D spatially anisotropic Heisenbcrg antifcrromagnet 
[AFzz with stronger exchange coupling along the c-axis 
than in the (a, b) planes]. Finally, the MOFFA, MOAFF 
and MOAAF phases are variations of the basic Kugel- 
Khomskii spin-orbital order || obtained by rotating the 
magnetic and orbital structure by tt/2. For the MOFFA 
phase at E z =0, the orbitals have a two-sublattice struc- 
ture in the (a, 6)-planes (x 2 — z 2 and y 2 — z 2 on the A- and 
B-sublattice, respectively). Along the c-axis strong an- 
tiferromagnetic spin-spin couplings are found, while the 
spin couplings in the (a, b) planes are ferromagnetic with 
a strength ~ rj. 

The anomaly occurs at the origin (E Z: rf) — (0,0) of 
the phase diagram: a 3D antiferromagnet (AFzz), a 2D 
antifcrromagnet (AFxx) and a quasi-lD A-type antifer- 
romagnet (MOFFA/MOAFF/MOAAF) become degen- 
erate! The emphasis on the 'uniaxial pressure' E z is mis- 
leading in the sense that the full scope of the problem 
is not visible directly from this phase diagram: at the 
origin of Fig. 1 an infinity of classical phases become de- 
generate. This is trivial to understand. In the absence of 
Hund's rule exchange, the Hamiltonian Eq. (||) becomes 
the full story. Assuming a 3D classical antiferromagnet, 

1) yields, 



Si ■ Sj = —1/4, and inserting this in Eq 



Heff = J L ^+ T t 



1 



(5) 



The orbital degrees of freedom are completely decou- 
pled and all 2 N orbital configurations have the same 
energy rf =0)! In addition, this infinity of dif- 
ferent 3D spin systems has the same energy as the 
MOFFA/MOAFF/MOAAF phases. It is actually so that 
at any finite temperature the 3D antiferromagnet be- 
comes stable because of the entropy associated with the 
decoupled orbital sector B. 

This 'gauge' degeneracy is clearly a pathology of the 
classical limit. We continued by studying the stability 
of the classical phase diagram with respect to Gaussian 
quantum fluctuations. As discussed in more detail in Ref. 
H this is a somewhat subtle affair. Intuitively, one could 
be tempted to think that the orbitals and spins can be 
excited independently. This is however not the case. The 
dynamical algebra of relevance to the problem is an so(4) 
algebra, and this implies that modes will occur which 
excite at the same time the spins and the orbitals: the 
spin-and-orbital waves (SOW)'s. Next to a (longitudinal) 
sector of pure orbital excitations, a 'transversal' sector 
is found corresponding with spin-excitations which are 
mixed with spin-and-orbital excitations, except for the 
acoustic modes at long wavelength which become pure 
spin-waves as imposed by the Goldstone theorem. 

We found that upon approaching the infinite critical 
point, the mass gap associated with the discrete symme- 
try in the orbital sector collapses. The (mixed) transverse 
modes give the dominating contribution to the renormal- 
ization of energy and magnetic order parameter. In the 
AFxx (AFzz) phase the lowest transverse mode softens 
along k = (n,0,k z ) [k — (fc x ,0, 0)], and equivalent lines 
in the Brillouin zone (BZ), regardless how one approaches 
the critical lines. Thus, these modes become dispersion- 
less along particular (soft-mode) lines in the BZ, where 
we find finite masses in the perpendicular directions, 
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with Aj = and ^ at the M point, and 
quantum fluctuations diverge logarithmically, (SS Z ) ~ 
/ d 3 k/uj(k) - / d 2 k/(A t + B t k 2 ) ~ In A. ( , if A, -> at 
the transition. We found that the quantum correction to 
the order parameter (S z ) becomes large, well before the 
critical point is reached. In Fig. 1 the lines are indicated 
where | (<5S'^) | = (S z ): in the area enclosed by the dashed 
and dotted lines classical order cannot exist, at least not 
in gaussian order. 

If the classical limit is as sick as explained in the previ- 
ous paragraphs, what is happening instead? A priori it is 
not easy to give an answer to this question. There are no 
'off the shelf methods to treat quantum spin problems 
characterized by classical frustration, and the situation is 
similar to what is found in, e.g. J1 — J2 — J3 problems [10|. 
A first possibility is quantum order-out-of-disorder [11|: 
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quantum fluctuations can stabilize a particular classical 
state over other classically degenerate states, if this par- 
ticular state is characterized by softer excitations than 
any of the other candidates. Khaliullin and Oudovenko 
Jl2| have suggested that this mechanism is operative in 
the present context, where the AFzz 3D anisotropic an- 
tiferromagnet is the one becoming stable. Their original 
argument was flawed because of the decoupling procedure 
they used, which violates the so(4) dynamical algebra 
constraints || . However, Khaliullin claims to have found 
an 'so(4) preserving' self-consistent decoupling procedure 



which does yield order-out-of-disorder |13|. Nevertheless, 
there is yet another possibility: valence-bond (VB) sin- 
glet (or spin-Peierls) order, which at the least appears in 
a more natural way in the present context than is the 
case in higher dimensional spin-only problems, because 
it is favored by the directional nature of the orbitals. 



The essence of a (resonating) valence bond [(R)VB] 
state is that one combines pairs of spins into singlets. 
In the short-range (R)VB states these singlets involve 
nearest-neighbor spin pairs. Subsequently, one particu- 
lar covering of the lattice with these 'spin-dimers' might 
be favored (VB or spin-Peierls state), or the ground state 
might become a coherent superposition of many of these 
coverings (RVB state). On a cubic lattice the difficulty 
is that although much energy is gained in the formation 
of the singlet pairs, the bonds between the singlets are 
treated poorly. Nevertheless, both in ID spin systems 
(Majumdar-Ghosh Q, AKLT-systems [pj) and in the 
large N limit of SU (N) magnets in 2D, ground states are 
found characterized by spin-Peierls/VB order fill . 



It is straightforward to understand that the interplay 
of orbital- and spin degrees of freedom tends to stabilize 
VB order. Since the orbital sector is governed by a dis- 
crete symmetry, the orbitals tend to condense in some 
classical orbital order. Different from the fully classi- 
cal phases, one now looks for orbital configurations op- 
timizing the energy of the spin VB configurations. The 
spin energy is optimized by having orbitals 3£ 2 — r 2 on 
the nearest-neighbor sites where the VB spin-pair lives, 
with £ directed along the bond. This choice maximizes 
the overlap between the wave functions, and thereby the 
binding energy of the singlet. At the same time, this 
choice of orbitals minimizes the unfavorable overlaps with 
spin pairs located in directions orthogonal to The net 
result is that VB states are much better variational so- 
lutions for the KK model, as compared to the standard 
Heisenberg spin systems. 
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FIG. 2. A variety of valence bond solids (see text). 

Adressing this systematically, we found that two fam- 
ilies of VB states are most stable: (i) The 'staggered' 
VB states like the PVBA and PVBIc states of Fig. 2. 
These states have in common that the overlap between 
neighboring VB pairs is minimized: the large lobes of the 
3£ 2 — r 2 wave functions of different pairs are never point- 
ing to each other, (ii) The 'columnar' VB states like the 
VBc (or VBa) state of Fig. 2. In the orbital sector, this 
is nothing else than the AFzz state of Fig. 1 (3z 2 — r 2 
orbitals on every site). Different from the AFzz state, the 
spin system living on this orbital backbone is condensed 
in a ID spin-Peierls state along the z-direction which is 
characterized by strong exchange couplings. The spins 
in the a(fe)-directions stay uncorrelated, due to the weak- 
ness of the respective exchange couplings as compared to 
the VB mass gap. 

The energies of these VB states and the classical states 
dressed up with quantum fluctuations are quite close to- 
gether. A key issue is if the true ground state is com- 
pressible (dressed classical state), or characterized by a 
dynamical mass-gap (VB states). This will most likely 
depend on subtleties beyond the reach of the relatively 
crude variational Ansatze presented here (l?]] . So the na- 
ture of the ground state of the Kugel-Khomskii problem 
for small Hund's-rule coupling is still an open problem. 



III. ELECTRON-PHONON COUPLING IN KCuF 3 

In the previous Section we discussed the orbital order 
as driven by the electron-electron interactions. However, 
one can think quite differently about the real systems: 
the deformations found in KCUF3 (or LaMnOa) could 
in principle be entirely caused by phonon-driven collec- 
tive Jahn- Teller effects. This subject has been intensely 
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studied in the past and is well understood. It starts out 
neglecting electron-electron interactions, and the focus is 
instead on the electron-phonon coupling. In case that the 
ions are characterized by a Jahn- Teller (orbital) degen- 
eracy, one can integrate out the (optical) phonons, and 
one finds effective Hamiltonians with phonon mediated 
interactions between the orbitals. In the specific case of 
e g degenerate ions in a cubic crystal, these look quite 
similar to the KK Hamiltonian, except that the spin de- 
pendent term is absent Jf8{ . Any orbital order resulting 
from this Hamiltonian is now accompanied by a lattice 
distortion of the same symmetry. 

The size of the quadrupolar deformation in the (a, b) 
plane of KCuF 3 is actually as large as 4 % of the lattice 
constant (a) . It is therefore often argued that the orbital 
order is clearly phonon-driven, and that the physics of the 
previous section is an irrelevancy. Although appealing at 
first sight, this argument is flawed: large displacements 
do not necessarily imply that phonons do all the work. 

The deformations of the lattice and the orbital degrees 
of freedom cannot be disentangled using general princi- 
ples: they constitute an irreducible subsector of the prob- 
lem. The issue is therefore a quantitative one, and in the 
absence of experimental guidance one would therefore 
like to address the issue with a quantitative electronic 
structure method. The LDA+U method is the method 
of choice. It is constructed to handle the physics of elec- 
tronic orbital ordering, keeping the accurate treatment of 
the electron-lattice interaction of LDA intact. According 
to LDA+U calculations the total energy gained by the 
deformation of the lattice is minute as compared to the 
energies involved in the electronic orbital ordering p9| . 
At the same time, the phonons are important on the 
macroscopic scale and they contribute to driving KC11F3 
away from the infinite-critical point of the phase diagram 
Fig. 1. 

We start out with the observation that according to 
LDA KCuF 3 would be an undistorted, cubic system: 
the energy increases if the distortion is switched on (see 
Fig. 3). The reason is that KCUF3 would be a band 
metal according to LDA (the usual Mott-gap problem) 
with a Fermi-surface which is not susceptible to a band 
Jahn- Teller instability. LDA+U yields a drastically dif- 
ferent picture Jh| . LDA can be looked at as unpolarized 
LDA+U, and by letting both the orbitals and the spins 
polarize an energy is gained of order of the band gap, 
i.e., of the order of 1 eV. The orbital- and spin polariza- 
tion is nearly complete and the situation is close to the 
strong coupling limit underlying the spin-orbital models 
of Section II. Also when the cubic lattice is kept fixed, 
the correct orbital and spin ordering (MOFFA of Fig. 1) 
is found, with spin-exchange constants which compare fa- 
vorably with experiment [^9). Because the orbital order 
has caused the electron density to become highly unsym- 
metric, the cubic lattice is unstable. Further energy can 
be gained by letting the lattice relax. The lattice distor- 



tion calculated in LDA+U (~ 3% of a) comes close to the 
actual distortion of KCuF 3 (~ 4 %). However, despite 
the fact that the distortion is large, the energy gained by 
the lattice relaxation is rather minute: ~ 50 meV (see 
Fig. 3)! Obviously, in the presence of the electronic or- 
bital order the cubic lattice becomes very soft with regard 
to the quadrupolar distortions and even a small electron- 
phonon coupling can cause large distortions. 
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FIG. 3. The dependence of the total energy of KCUF3 on 
the quadrupolar lattice distortion according to LSDA and 
LDA+U band structure calculations (after Ref. 

Although the energy gained in the deformation of the 
lattice is rather small, the electron-phonon coupling is 
quite effective in keeping KCuF 3 away from the physics 
associated with the origin of the phase diagram (Fig. 1). 
Since the ferromagnetic interactions in the (a, b) plane of 
KCuF 3 are quite small ( J a b = — 0.2 meV, as compared to 
the 'ID' exchange J c = 17.5 meV [M), one might argue 
that the effective Hund's rule coupling Jrj as of relevance 
to the low energy theory is quite small. Although this still 
needs further study, it might well be that in the absence 
of the electron-phonon coupling KCuF 3 would be close 
to the origin of Fig. 1. However, the electron-phonon 
coupling can be looked at as yet another axis emerging 
from the origin. In principle, the electron-phonon cou- 
pling introduces two scales: (i) a retardation scale, which 
is governed by the ratio of the phonon frequency and the 
electronic scale set by J ~ 20 meV. Since J is relatively 
small, KCuF 3 is close to the anti-adiabatic limit where 
the lattice follows the electronic fluctuations, (ii) in the 
anti-adiabatic limit the phonons are high energy modes 
which can be integrated out, causing the effective orbital- 
orbital couplings we earlier referred to. These couplings 
destroy the cancellations leading to Eq. (|)|), thereby driv- 
ing the system away from the point of classical degener- 
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acy. The typical scale for the phonon induced effective 
orbital interactions is at most of the order of the LDA+U 
lattice relaxation energy. However, as the latter (~ 50 
meV) is quite a bit larger than J, the effective interac- 
tion will likely be able to put KCuF 3 well outside the 
'dangerous' region near the origin of the phase diagram. 

In summary, although further work is needed it might 
be that phonons are to a large extent responsible for the 
stability of KCuFa's classical ground state. In any case, 
one cannot rely on the sheer size of the lattice deforma- 
tions to resolve this issue! 



IV. HOW ABOUT THE MANGANITES ? 

Given the discussion so far, the search for interest- 
ing quantum effects in orbital degenerate Mott-insulators 
should not be regarded as hopeless. Unfortunately, the 
insulating parent compounds of the CMR manganites, 
such as LaMnC>3, are not candidates for this kind of 
physics. The reason is not necessarily phonons: also 
in the manganites the 'Jahn- Teller' lattice distortions 
are sizable, but this does not necessarily imply that the 
phonons are dominating. Two of us derived a Kugel- 
Khomskii-type model of relevance to this regime, and we 
did find a dynamical frustration of e 9 -superexchange at 
Jh — pi] ]. However, the system is driven away from 
this point by two effects: (i) the manganites are in the 
Hund's rule dominated regime, with a large splitting be- 
tween the lowest energy high-spin state at U — 5Jh (with 
Jh = 0.69 cV p3|), and the low-spin states at energies 
~ U ; (ii) the additional ^g-superexchange between the 
5 = 3/2 cores favours an antiferromagnetic order in all 
three spatial directions. The net outcome is that the fer- 
romagnetic interaction between the total S = 2 spins in 
the (a, b) planes is of order of the c-axis exchange, sig- 
nalling that the manganites are in the Hund's rule stabi- 
lized regime of the phase diagram. 

The mysteries of the manganites relate to what hap- 
pens when quantum-mechanical holes are added to the 
orbital/spin ordered insulator. This is undoubtedly a 
problem with its own characteristics, which cannot be 
reduced to a variation on the far simpler problems en- 
countered in the insulators. Nevertheless, we do believe 
that the study of the insulating limit might be of some 
help in better appreciating what is going on in the doped 
systems. It is tempting to think about orbital degrees of 
freedom as being spins in disguise. This is not quite the 
case. Orbitals are far less quantum-mechanical - they are 
more like Ising spins than Heisenbcrg spins. Secondly, or- 
bitals carry this unfamiliar property that depending on 
their specific orientation in internal space, overlaps in- 
crease in particular real space directions, while they di- 
minish in orthogonal directions. Our valence-bond con- 
structions illustrate this peculiar phenomenon in the case 
of spins, but the same logic is at work when the hole is 



delocalizing. This intimate connection between internal 
symmetry and the directionality of derealization causes 
the dynamical frustration which has been highlighted in 
this communication. This motive seems also at work 
in the doped system, witness the many near degener- 
ate states found both in mean-field calculations |p2 24 
and in experiment ]25[ ]. Further work is needed on this 
fascinating problem. 
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